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Abstract 

The evolution of the logarithmic scaled energy spectra with energy scale, in the framework of the modified 
leading logarithmic approximation (MLLA), is investigated. The behaviour of the higher order moments of 
HERA deep inelastic scattering data are compared with the analytic predictions of MLLA. 
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1 Introduction 



The properties of the hadronic final state cannot be calculated dhectly in Quantum Chroniodynamics (QCD). 
It is necessary to model the transition from the perturbatively calculable partonic final state to hadrons. One 
approach that has been successful in describing the general features of the inclusive energy spectra, in both e+e" 
annihilation and deep inelastic scattering (DIS) experiments [0, is the modified leading log approximation 
(MLLA). This uses the local parton hadron duality (LPHD) approach |^ to relate the MLLA partonic predictions 
to the hadronic observables. 

The hadronic final state reflects the underlying partonic behaviour. At small values of momentum fraction 
carried by the outgoing parton from the original hard scatter, the parton shower development is dominated by 
gluon bremsstrahlung. The partonic properties can be calculated, on average, in the framework of the MLLA. 
These MLLA calculations have two free parameters: a running strong coupling, defined by a QCD scale A, and 
an energy cut-off, Qoi below which the parton evolution is truncated. 

LPHD assumes that a simple normalisation factor, which does not explicitly include either the effect of 
particle masses or resonance decays, can be used to relate the hadronic distributions to the partonic ones. At 
large enough energies, away from where the parton evolution is truncated, this 'hadronisation' factor should be 
approximately independent of the energy scale ^ at which the parton spectra are calculated. 

In our previous study Q], various approaches were taken to calculate the single particle spectra and their 
moments within the MLLA framework and compared with e^e~ data. In this paper the theoretical results are 
compared with deep inelastic scattering data in the current region of the Breit frame . 

The event kinematics of DIS are determined by the negative square of the four-momentum of the virtual 
exchanged boson, = —q^, and the Bjorken scaling variable, x — /2P-q, where P is the four-momentum of 
the proton. In the Quark Parton Model (QPM), the interacting quark from the proton carries four- momentum 
xP. 

A natural frame in which to study the dynamics of the hadronic final state in DIS is the Breit frame 1^. In 
this frame, the exchanged virtual boson is completely space-like and has a four-momentum q = (0,0, 0,— Q = 
—2xP^^^^*), where p^'"'^** is the momentum of the proton in the Breit frame. The particles produced in the 
interaction can be assigned to one of two regions, known as the current and target regions. The Breit frame 
allows the maximal separation of the incoming and outgoing partons in the QPM. 

The current region in the Breit frame of DIS is analogous to a single hemisphere of e+e^ annihilation. In 
the current region, the outgoing struck quark from the proton has a momentum —Q/2. The scaled momentum 
spectra of the particles in the current region, expressed in terms of — ln(Q/2p^'''^**), are expected to have 
a dependence on Q similar to that observed in e+e~ annihilation [Q at energy .^/s^ = Q. It should be noted 
though, that the current region of DIS will have a different primary quark flavour composition than e+e~ 
annihilation, in particular for heavy quarks. 

This paper extends the original ZEUS analysis by investigating the effect of introducing an effective mass 
term into the MLLA predictions (see section |[) In contrast to the previous analysis of c+c^ data Q this 
effective mass was treated as an additional free parameter of the theoretical models and was fitted to the data. 



2 Single Particle Spectra 

Given a high-energy parton that fragments via secondary partons into a jet of hadrons, the MLLA evolution 
equation allows the secondary parton spectra for the logarithmic scaled energy, ^, to be calculated 0. The 
variable ^ is defined as ln(_Eo/-E) = ln(l/a;), where Eo{= Q/2 for DIS) is the energy of the original parton and 
E is the energy of the secondary parton. The cut-off, Qq, limits the parton energy io E > kx > Qo, where kx 
is the transverse energy of the decay products in the jet evolution. In order to reconstruct the ^ distributions 
an inverse Mellin transformation is performed. 

It is convenient to investigate the MLLA spectra in terms of moments. The analytic form of the cumulant 
moments has been calculated in for the first four moments. This allows the normalised moment, (^™), to be 
calculated and hence the dispersion (cr), skewness (s) and the kurtosis (fc) to be constructed, see for example 



(1) 



a' - if2 = r>-(er (2) 

k = Ki/a = J . (4j 

The so-called limiting spectrum [|j is the case when Qq = A. The ^ distribution, F), is dependent on 
one parameter Y, which is given by ln(i?o/A), where A is now the only free parameter. For this limiting case, 
Fong and Webber |^ have also calculated the behaviour of the moments of the ^ spectra with energy scale Q. 
They point out that the spectra can be represented close to the maximum of the distribution by a distorted 
Gaussian of the form: 



D{^,Y) cx exp 



(5) 



where 5 = {^ — l)/a and / is the mean. The differences between the parton spectra for quark (and gluon) jets are 
Next-to-MLLA effects and have been included in these investigations according to the approximate prescription 
of Dokshitzer et al. ]Tot . 

The MLLA limiting spectra predicts that the i?o-dependence of ^max should have the form 

?n.a. = rf0.5+./g-g+f4V\ (6) 



Y Y \Y^ 

where C is expected to be w 0.3 [ pl| . 

The mean and peak position, ^max, are related, for the skew Gaussian used here, by 

emax=(0-- + -(l+-^')^- (7) 
(7 S 

Thus the maximum is a function of several variables. 

In the MLLA approach discussed above, the partons are assumed massless so the scaled energy and mo- 
mentum spectra are identical. Experimentally the scaled momentum distribution is usually measured and as 
the observed hadrons are massive the equivalence of the two spectra no longer holds. In JT2[ the assumption is 
made that the cut-off Qq can be related to an effective mass, rricff, of the hadrons. This allows the logarithmic 
scaled momentum distribution, ^p, to be written as 

where 

^ = log ■ 



and the energy of a hadron with a momentum ph is Eh — \/pfi + "^cff- "^^^ limiting momentum spectra based 
on massless partons or massive partons with the uics being set to Qo = A, will be referred to as MLLA-0 and 
MLLA-M spectra, respectively. 



3 Comparison with DIS Data 



Figure ^ shows the differential £,p distributions, with statistical errors only, for charged particles in the current 
fragmentation region of the Breit frame for nine (Q) values from the ZEUS experiment These distribu- 
tions are approximately Gaussian in shape with the mean charged multiplicity given by the integral of the 
distributions. As Q increases, the multiplicity increases and, in addition, the peak of the distribution moves 
to larger values of ^p. The data were fitted to the distorted Gaussian of eqn. ^ over a range of ±1.5 units (for 
(Q) < 14.5 GeV) or ±2 units (for (Q) > 14.5 GeV) in around the mean. This is the range used by the ZEUS 
collaboration in their analysis of the data. The smooth curves in Fig. |^ result from the fit to the data; they 
represent the data well. The fitted parameters were insensitive to the range of the fit. The values are not 
quoted here because the measurements are dominated by systematic errors. 

The ZEUS data, solid points, are compared with various MLLA curves in Figure |[ For clarity both the data 
and theory are normalised to a maximum value of unity. The theoretical curves (solid line MLLA-M, dotted line 
MLLA-0) are calculated for A = 280 MeV(= TOcff)- Both the MLLA-0 and MLLA-M describe well the general 
trend of the data. In particular, MLLA-M describes well the central region of data for (Q) > 10.4 GeV. There 
is a trend however, for the MLLA-M distribution to be more positively skew than the data; this is particularly 
pronounced at low energies. In contrast, MLLA-0 is consistently less skew than the data with the skewness 
being negative at low Q. MLLA-M and MLLA-0 converge to similar distributions as Q increases. At low Q 
the data are peaked at lower values of than the theoretical predictions while at high Q the predictions peak 
lower than the data. 

The histogram is a fit of MLLA, in the theoretically allowed region, over the range of data above the 
horizontal dashed line (80% of the data as measured from the maximum of the distributions), taking both A 
and TOcff to be free parameters. The normalisation of the distribution in each Q bin was allowed to vary. The 
fit gives A — 280 ± 2 MeV and rrics — 230 ± 6 MeV. Again, within the errors, the parameters were insensitive 
to the region of the fit to the data. The fitted MLLA gives a distribution intermediate between MLLA-M and 
MLLA-0. The data are well fitted for (Q) > 10.4 GeV. Otherwise, the fit closely follows MLLA-M and shows a 
similar disagreement with the data. 

To study further the differences between the MLLA predictions and the data, fits of a skew Gaussian have 
been made to the MLLA predictions, in a manner similar to the fits made to the data, in order to extract the 
moments. Figure ^ shows the first four moments of the spectra as a function of Q. For the ZEUS data (solid 
points), the mean and a increase as Q increases, whilst the skewness and kurtosis fall. Also shown are the 
results from the fits to the MLLA-M and MLLA-0 spectra, the upper and lower shaded areas respectively. For 
both MLLA-0 and MLLA-M the range fitted lay between 30% of the spectrum (lower limit of band) to 80% of 
the spectrum (upper limit of band). The upper limit of the fit in ^p was determined by the limit of the MLLA 
spectrum if this were less than the percentage range. In all cases the skew Gaussian describes the theory to 
better than 0.5%. 

There are two points to be noted when making these fits: firstly the fitted parameters of the theory depend 
on the range of spectrum fitted. Secondly, fitting over a range greater than ±la may mean that the skewness 
and kurtosis are not those of the distribution because the approximations made in deriving the expression 
become invalid. Consequently the the fitted parameters should be considered as simply a parameterization 
of the distributions rather than accurate measures of the skewness and kurtosis. Equally because the a and 
kurtosis are correlated, the measured a is fit-range dependent. 

For cr, skewness and kurtosis, Fig. ^, the DIS data are bracketed by the MLLA-M and MLLA-0 predictions 
as would be expected from the plots of Figure |[ The skewness predicted by MLLA-M is more positive than 
the data while MLLA-0 underestimates the skewness; for the highest Q both models approach the data. There 
is evidence of a systematic difference in the Q-dependence of the mean between the models and the data for 
the energy range studied here. 

The dotted line shows the moments obtained from a fit to the massive MLLA spectrum obtained using the 
fitted A = 280 MeV and TOoh = 230 MeV. It is evident that a reasonable description of all moments other than 
the mean is obtained; the mean has an energy dependence that differs from that of the data. 

The energy dependence of the ^p distribution has been explored further by studying the dependence of the 
peak position on Q. This has the advantage, in comparison with the moments, that the maximum is relatively 
insensitive to the function used to parameterize the shape of the S,p distribution. The Q-dependence of ^max is 
shown in Figure ^. The shaded bands show the maxima derived from MLLA-0 and MLLA-M. For both cases 
the range of fit is from 30% to 80%. As before the dotted line uses the MLLA parameters A = 280 MeV and 



rrieff = 230 MeV from the fit. This plot confirms that the energy dependence of the DIS data differs from that 
of any of the MLLA predictions. 

A good fit of eqn. ^ to the ^max from the experimental data can be achieved with A = 285 MeV (with d set 
to zero), which is consistent with the A value from the fit to the spectrum. However the value of C is a factor 
two larger than expected and indicates again an inconsistency of the energy dependence of the data with the 
expectations of the MLLA. 

These results are consistent with our previous analysis of e+e~ data [Q. 

4 Conclusions 

The HERA DIS data in the current fragmentation region of the Breit frame have been compared to the MLLA 
predictions as a function of Q. In general, the data are well described within the LPHD-MLLA approach, with 
the predictions of MLLA-0 and MLLA-M bracketing the data. With increasing Q, the various predictions of the 
MLLA converge on the data. A fit has been made to the data using the massive MLLA approach and values of 
A and rrics are found that give a reasonable description of the moments of the ^ distribution. 

In detail discrepancies are found between the theoretical predictions and the data. In particular, a consistent 
description of the evolution of the peak position of the ^ distribution as a function of scale Q is not possible; 
no energy-independent combination of A and mcs is found that can describe the data. 
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Figure 1: The charged particle distributions of 1/Ndn/d^p in the current fragmentation region as a function of 
for different Q bins. Only statistical errors are shown. The full line is the skew Gaussian fit. 
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Figure 2: The charged particle distributions of 1/Ndn/d^p in the current fragmentation region as a function 
of for different Q bins. Only statistical errors are shown. The full and dotted lines are the MLLA-M and 
MLLA-0 predictions with A — 280 MeV respectively. The histogram is fit of the MLLA to the data with 
A = 280 ± 3 MeV and m^g = 230 ± 8 MeV. 
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Figure 3: The moments of the distribution in the current fragmentation region as a function of Q, from fits to 
ZEUS data, are shown as solid points. The errors are those returned from the fit. The upper and lower shaded 
areas are from fits to the MLLA-M and MLLA-0 spectra, respectively. The dotted line shows the moments 
obtained from a fit to the MLLA spectrum obtained using the fitted A = 280 MeV, nics = 230 MeV. 
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Figure 4: The peak position of the distribution in the current fragmentation region as a function of Q. The 
errors shown are those returned from the fit. The upper and lower shaded areas are from fits to the MLLA-M 
and MLLA-0 spectra, respectively. The dotted line shows the ^max obtained from a fit to the MLLA spectrum 
obtained using the fitted A = 280 MeV, Wcff = 230 MeV. 



